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Electron counting of a single porphyrin molecule between two electrodes shows a crossover from 
sub- to super-Poissonian statistics as the bias voltage is scanned. This is attributed to the si- 
multaneous activation of states with electron transfer rates spanning several orders of magnitude. 
Time-series analysis of consecutive single electron transfer events reveals fast and slow transport 
channels, which are not resolved by the average current alone. 



Real-time measurements of single-electron trans- 
fers through nanosystems have been reported 
recentlj^i^ii^i^i^i^. Although electron counting has 
not yet been observed in single molecule junctions, new 
techniques based on carbon nanotubcs have been pro- 
posed for its possible realization^. Theoretical efforts are 
required to fully connect the information of the electron 
transfer statistics to the molecular properties. This 
will offer new opportunities to study many-body effect 
to characterize the bonding of molecules to external 
electrodes with implications to molecular electronics^i^. 
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counting 
provides 

detailed information about the probability P{k, t) of 
transferring a net-number k = kin — kout of electrons 
between electrode and molecule during a time interval 
t. The cumulants of P{k,t) are directly related to 
important properties of the junction. For example, 
the first order cumulant Ci(i) = k/t gives the average 
current I{t) — eCi{t). The shot noise is related to the 
second cumulant S{t) = 2e'^C2it) = 2e'^(k^ - k'^)/t 
and is commonly represented by the Fano fac- 
tor F{t) = C2{t)/Ci{t). The third cumulant 

Cj,{t) = {k — kY/t measures the skewness of the 
distribution. 

The goal of the present work is to relate the higher sta- 
tistical cumulants to the intrinsic properties of the metal- 
molecule interface which cannot be accessed by average 
current measurements. We simulate a magnesium por- 
phine (MgP) molecule coupled to two metallic electrodes 
as shown in Fig.[T] and find that the Fano factor shows a 
crossover from sub- {F < 1) to super-Poissonian (F > 1) 



shot noise when a large number of transport channels are 
opened by the bias. 

Mechanisms like asymmetric coupling strengths of the 
two clectrode o^^'^^1^^'^'^ , Coulomb charging energy in 
multi-level systems^^iSl and potential fluctuations in res- 
onant quantum wells^^ were predicted to induce super- 
Poissonian transfer statistics. In case the two electrodes 
are coupled to the molecule with equal strength, we iden- 
tify an additional mechanism in the MgP-metal bonding 
and attribute it to the presence of couplings between the 
electronic many-body states of the molecule spanning dif- 
ferent orders of magnitude. This is verified using a simple 
model. 

The presence of these overlapping effects makes it dif- 
ficult to analyze specific electronic states and couplings 
using the statistical cumulants. We further examine the 
probabilities of consecutive electron transfers which can 
be obtained from a time-series of single electron transfer 
events^^ and calculate a decay time distribution of the 
MgP molecule junction. We demonstrate that slow and 
fast transport channels can be unambiguously identified 
as peaks in the decay time distribution. 

The junction Hamiltonian is decomposed as H = Hq + 
Hs + Hd + Ht, where Hq represents the isolated MgP 
molecule, Hs(^]j^ the source (drain) electrode, and Ht 
is the electrode-molecule interaction. We assume that 
the density of states of the electrode is constant on the 
relevant energy scale and that the current is small. 

Describing the electrodes as free electron reservoirs and 
neglecting coherence elements, the Pauli rate equation 
provides a convenient framework to include many-body 
effect s— i-S and electronic structure calculations'^"'^. The 
Pauli equation for the population pl of state \i), eigen- 
state of Hq, with energy Ei and TV electrons is then given 

by 
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The molecule can be excited by electron transfers from 
the electrodes. Only transitions between states with a 
unit charge difference are allowed. The transition rates 
Wi-,j involve contributions from source {a = S) and 
drain {a — D) electrode and may be decomposed as 
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WlZ'j = Wi::f' + Wl^f where W^ 
is the rate of electron transfer into 



the 



molecule and 
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,1 — fa{Eij)) accounts for the reverse pro- 
jj — Ei — Ej is an energy difference between 
charge states and fa{E) = (1 + exp[f3{E — ^1^)])^^ is 
the Fermi function for a free electron gas with thermal 
energy 1/P and chemical potential /i^. Each pair W~^^ 
of transition rates defines a transport channel from the 
source to the drain which can be opened and closed by 
the chemical potentials. The transition rate between a 
pair of states depends on the electrode-molecule inter- 



action through the couplings Tlj = J2s 1"^^ \'^\{j\cs\i)\'^ 
where the index s runs over the molecular orbitals and 
Ts is a coupling strength parameter weighted by the 
overlap {j\cs\i) between many-body states. 

Ground and excited states of the neutral and the 
charged MgP are calculated at the Hartree Fock (HF) 
and configuration interaction singles (CIS) level using a 
6-31G basis set'^^. We fix the molecular geometry of the 
minimum energy configuration of the neutral MgP. Fig. [2] 
shows the nine lowest electronic energies of the neutral 
N, anionic A^ -I- 1 and cationic N — 1 MgP charge states; 
the neutral molecule's ground state energy is taken to be 
zero. The connecting lines are a guide for the eye. These 
states are connected by a network of couplings r^"\ For 
states |i) and \j) the overlap factor {j\cs\i) is calculated 
using the CIS expansion coefficients, the molecular or- 
bital coefficients, and the atomic overlap matrix; the pro- 
cedure is described elsewhere^. The overlap factors for 
MgP show that multiple orbitals contribute to the elec- 
tron transfer rate indicating that a many-electron picture 
is necessary to describe the transport. These are shown 
in the supplementary information section. 

The Ts factors are assumed to take the form of 
a kinetic energy integral between molecular orbital \s) 
and a spherically symmetric probe orbital \a) centered 
about a point r^ representing the lead-molecule con- 
tact. The probe orbitals used in our calculations are 
centered at the tips of the gray paraboloids shown in 
Fig.[T] A useful representation of the coupling is obtained 
by assuming the probes have perfect spatial resolution, 
i.e., (r\a) = S{r — Tq). In this case Ts' ex V'^(j)s{ra) 
where 4)s{fa) is the orbital wavefunction evaluated at 
the probe's center such that the coupling Ff- can be 
readily visualized in 3D space. In Fig.[3l we show the 
Fii, Fig, Fig, F21, F23 and F28 couplings between the 
neutral and anionic charge state plotted with respect to 
a delta-function probe's central position. The couplings 
are highly sensitive to the electrode-molecule orientation. 



The steady state is calculated by setting pi ' = in 
Eq. 11]). From the resulting homogeneous system of equa- 
tions an inhomogeneous form can be derived using the 
additional constraint "Y^^Pi — 1. For an irreducible rate 
matrix, the inhomogeneous system can be solved alge- 
braically for p. 

We use a generating function technique to cal- 
culate the electron transfer statistics through the 
junction^^'ii^'i^^'i^^^^^^^^. The asymp- 
totic cumulants, describing an infinite time counting 




FIG. 1: The geometrical configurations (a),(b) and (c) of the 
source (top) and drain (bottom) electrode with respect to the 
MgP molecule. 




• neutral (N) 
■ anion (N+1) 

♦ cation (N-1) 



2 4 6 8 

electronic quantum number 

FIG. 2: The eigenenergies of the electronic states for the neu- 
tral, anionic and cationic charge states of the MgP molecule. 



measurement, can be calculated from the eigenvalues of 
the generating function propagator^^-^^. For details see 
supplementary material. The statistics is generated by 
counting the net-transfer (in minus out) at the source. 
Since [3V > 1 in our calculations, electron transfers in 
opposite direction of the bias are insignificant and the 
net-process is practically reduced to transfers with the 
bias. 

We employ the electrode configurations (a),(b) and (c) 
shown in Fig.[TJ The probe orbitals are centered at the 
tips of the gray paraboloids. The upper electrode acts 
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FIG. 3: The spatial profiles of the couplings F 
electronic states {i,N) and {j,N + 1). 
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FIG. 4: The net-current through the MgP junction as a func- 
tion of voltage for electrode geometries (a), (b) and (c). The 
corresponding Fano factor and the skewness of the electron 
transfer statistics at the source electrode are quantified by 
the right scale. 



as the source, the bottom is the drain. The chemical 
potential of the source is fixed at ns — Ep ~ 2.{)eV and 
the drain chemical potential /i£) — Ep — V is decreased. 
The voltage is varied in the energy range of the employed 
states between and 4.75^. In this range we can safely 
limit our discussion to the neutral and anion states since 
the cationic states are inaccessible. The temperature is 
set to T = 50K. The current is given in nano Ampere 
(nA) and time in nano seconds (ns). 

In configuration (a) both electrodes are symmetri- 
cally located at carbon ring. The corresponding current- 
voltage curve (Fig. [4^) shows three steps when charge ex- 
citations in the molecule are energetically allowed. Since 
the equilibrium Fermi energy is set to Ep = 2.QeV, the 
first transition, between the neutral ground state N,i = 1 
and the lowest anionic state iV -I- l,i = 1, is shifted by 
2.QeV as well. Below a bias of 3.7eT^ only these two states 
contribute to the current. Higher bias causes a cascade of 
excitations and additional states become occupied. The 
population dynamics is shown in the supplementary ma- 
terial. 

Interesting features can be observed in the higher cu- 
mulants of the electron transfer statistics through the 
contact. Fig.H] shows the Fano factor and skewness 
C3/C2 on the right scale. For small bias, the current is 
small and the single electron transfers are uncorrelated 
leading to a Poissonian process with Fano factor F — 1. 
It drops to a sub-Poissonian regime -F < 1 when only 
the Fii coupling is active. However, it becomes super- 
Poissonian _F > 1 at higher bias. The skewness is even 



more sensitive and reaches negative values in the super- 
Poissonian regime &iV > i.SeV. 

It has been shown2ii2^iS&i22 that a strongly asymmetric 
coupling of the system to the source and drain electrode 
can cause super-Poissonian electron transfer. This can 
be ruled out since the junction is symmetrically coupled 
to drain and source. We find that the sub-to super- 
Poissonian crossover can be attributed to two mecha- 
nisms. One is well known due to Coulomb charging en- 
ergy in multi-level system o '^ . 

The second effect can be attributed to large or small 
ratios of active excitation rates W^'^jW^'^i from a given 
(N,i) state to two (N-|-l,j),(N-|-l,l). Because (3 is much 
smaller than the energy gaps between the electronic 
states, the Fermi function approximates a heavy-side step 
function and the ratio of the rates is directly related to 
the the corresponding couplings Ti^j /Ti^i. 

This generic effect is rationalized in the supplemen- 
tary material using a simple model. For ratios of the 
couplings of the same electrode larger than 5 or smaller 
than 1/5 we observe super-Poissonian electron transfer 
statistics. Since it does not depend on an asymmetric 
coupling between source and drain it can be considered 
as an intrinsic characteristic of the molecule-electrode in- 
terface. For example in Fig. [3] we see that F^g coupling 

is weak compared to F^"' at the edge of MgP and we find 
a ratio of Fj^" /Fj^g « 10 at the position of the tip. 

We next move both electrodes to the center of MgP 
as shown in Fig.[T] (b). The current, Fano factor and 
skewness are shown in Fig.|3}D. Again for a bias below 
3.7y only the M^fl^i channel is active and the electron 
transfer is sub-Poissonian as indicated by the Fano factor. 
Once the additional channels are opened by the bias the 
Fano factor increases but does not exceed unity since the 
coupling ratio of the most significant channels is within 
the limits for sub-Poissonian statistics as predicted by 
our model calculations. In the center of the MgP, F^g 
is strong as shown in Fig Fig. [3] and we have a ratio of 

1^11/^18 ~ 1/4 ^^ ^^^ ^iP- Here the increase of the Fano 
factor is due to the Coulomb charging and this electrode 
configuration provides a good quantitative estimate for 
it. 

Let us move the drain back to its position in configu- 
ration (a) while the source is kept in the center as shown 
in Fig.[T] (c). Besides the Coulomb charging and extreme 
ratios of the couplings, we now also have asymmetric cou- 
pling to drain and source and three possible sources for 
super-Poissonian transfer are present. However the spa- 
tial asymmetry can be estimated from the average cur- 
rent signal of the two symmetric configurations. The 
Fano factor in Fig. St shows the crossover but is with a 
maximal value of -F = 1.26 smaller than in configuration 
(a). Qualitatively, current signal, Fano factor and skew- 
ness resemble an average over the corresponding signals 
of (a) and (b) . The asymmetry effect is therefore weaker 
in configuration (c) than the other two. 

We next examine the elementary probabilities of con- 




then given by 
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FIG. 5: The decay time distribution Gs^oi'y) for the elec- 
trode geometries (a)-(c). A finite width of the peaks was 
introduced by modulating the signal with a frequency of 
Lo — 7.510*l/ns. A logarithmic 7-axis is used for the left 
side of the middle panel and the bottom panel. 



secutive electron transfer events^^^^^^. An electron is 
detected when it enters the junction through the source 
at time to = and leaves through the drain at time t. No 
other transfer events take place in between. We denote 
this path by 5 — > D. The two corresponding electron 
transfer operators are M^('^'+) and VJ/(^'~) respectively. 
Their matrix elements are given by the transfer elements 
of the Pauli equation ([T]). The elementary probability of 
this two-time measurement is then given by2^ 



P.' 



S^D 



(i,to) = (T4^(^-)t/(t,to)M/(^'+^ 



(2) 



where U{t,to) is the propagator of the system in ab- 
sence of transfer events at the electrodes within the time 



interval U(t,to) 



:,2^ci = S 
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The opera- 



tor M^("''''"3) contains diagonal elements of the Pauli 
rate matrix calculated by VF^"'''*"^) = J2i(^^)[Wl 
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In order to measure Ps^oit, to) one has to detect sin- 
gle directionally-resolved electron transfers between the 
electrodes and the system and record a sufficiently long 
time-series of transfer events. Then a histogram of the 
number of consecutive transfer events S ~* D as function 
of increasing time intervals t — tg can be generated. By 
normalizing the histogram by the total number of events 
in the time series one obtains Ps^_D(t, to). 

We can calculate the two point-probability in Eq.® 
analytically. We consider the electronic states |1) to \m) 
which belong to charge state N and the electronic states 
\m + 1) to \n) that belong to charge state N + 1. The 
elementary probability of the S ~>^ D transfer path is 



Ps^Dit,to)^ J2 e-'^^'-'«^AiBi. (3) 
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The decay rates 7; are given by 7; — Yl\ 

fD{Eu)) +T^i\l - fs{Eu)) only. The coefficients Au 
are determined by the source and drain couplings and the 
population of the states: Ai = X^fcLi ^ik (^ ~ foiEik)) 

^i = J2T=i ^ij Js{Eij)pj. The populations are given by 
the steady state of the system. Surprisingly, each decay 
rate includes coupling elements to a single state of charge 
state N only. 

In order to analyze the different time scales in 
Ps^oit, to) due to the exponents, we use a one sided in- 
verse Laplace transformation to calculate the decay time 
distribution Gs^oil) = /o°° dt e('^+^")*Ps^£>(t). This 
gives 



Gs^nh)^ Y. 



AiBi 



l=m+l 



7; - 7 + ZCJ 
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An inverse Laplace transformation over a multiple expo- 
nential decays of the form ([3]) results in a sum over delta- 
functions. The parameter uj was introduced for assure a 
finite width of the peaks in \Gs^D(n)\- The magnitudes 
Gs^D (7) I are shown in Fig.[n]for the three electrode con- 
figurations. We use a A.75V bias where a large number 
of channels is activated. Each peak in Fig.[5]is due to the 
different time scale occurring in Ps^oit)- The positions 
of the peaks are given by the sum over the magnitude of 
active couplings 7^ to a specific electronic state i of the 
anionic charge state TV -t- 1 and are identified in Fig. [5] 
For all three configurations, only the 71 peak related to 
the Fii coupling is visible at a bias of 3V. Its position 
is shifted with increasing bias as soon as additional cou- 
plings Tji to state (1, iV) are activated. The spectrum of 
electrode configuration (a) shows dominant peaks for the 
first three electronic states of the anion. The 73 peak is 
weak as one could expect by comparing with the spatial 
profile of F18 in Fig|31 On the contrary a strong 73 peak 
is found in configuration (b) where the electrodes are in 
the center of the MgP molecule. Generating a spectrum 
for different bias and Fermi level by recording the po- 
sition and the magnitude of the peaks enables one to 
resolve individual transport channels. The average cur- 
rent solely depends on the sum over the active couplings 
and one cannot tell them apart. Rastering the molecule 
with the electrodes would then allow to measure the spa- 
tial profiles of the couplings shown in FiglJl Note that if 
one applies a high bias to the source electrode, its Fermi 
function would simplify to fs{Eii) — 1 for all eigenstates 
and the decay rates 7; would depend on the properties of 
the coupling to the drain only. In this case the molecule 
could simply sit on a surface acting as electrode and a 
tip acting as drain could measure the spatial profile of 
the couplings. Of interest is also the decay time distribu- 
tion of transfers in opposite direction of the bias since the 



sum in 7; over the neutral electronic states would then 
be replaced by a sum over the anionic states. Combining 
both spectra would then provide detailed information on 
each coupling element F^ . Bi-directional electron trans- 
fer has been observed in quantum dots-. The practical 
limitations of the proposed method are set by the in- 
verse Laplace transformation which has to be perfomed 
numerically for experimental data. 

In summary, the second and the third cumulants pro- 
vide detailed information and can be used to probe ra- 
tios of the couphngs between electronic states due to the 
molecule-electrode bonding. This effect is explained us- 
ing a simple model. The decay time distribution of two 
point electron transfer probabilities can resolve active 



transport channels and the magnitude of the couplings 
between the electronic states. The prestent methods are 
general and applicable to different kinds molecules as well 
as quantum dots. 

Supporting Information Available: Additional 
theoretical and computational detail are presented. This 
material is available free of charge via the Internet at 
http://pubs.acs.org, 
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